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Abstract 


17 Global atmospheric energy balance is one of the fundamental processes for the earth’s 

18 climate system. This study uses currently available satellite data sets of radiative energy at the 

19 top of atmosphere (TOA) and surface and latent and sensible heat over oceans for the year 2000 

20 to assess the global annual energy budget. Over land, surface radiation data are used to constrain 

21 assimilated results and to force the radiation, turbulent heat, and heat storage into balance due to 

22 a lack of observation-based turbulent heat flux estimations. 

23 Global annual means of the TOA net radiation obtained from both direct measurements 

24 and calculations are close to zero. The net radiative energy fluxes into the surface and the 

25 surface latent heat transported into the atmosphere are about 113 and 86 W/m~, respectively. 

26 The estimated atmospheric and surface heat imbalances are about - 8 ~ 9 W/m~, values that are 

27 within the uncertainties of surface radiation and sea surface turbulent flux estimates and likely 

28 systematic biases in the analyzed observations. The potential significant additional absorption of 

29 solar radiation within the atmosphere suggested by previous studies does not appear to be 

30 required to balance the energy budget the spurious heat imbalances in the current data are much 

31 smaller (about half) than those obtained previously and debated at about a decade ago. Progress 

32 in surface radiation and oceanic turbulent heat flux estimations from satellite measurements 

33 significantly reduces the bias errors in the observed global energy budgets of the climate system. 


34 
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Introduction 


36 Global atmospheric energy and heat balance is one of the fundamental physical processes 

37 of the earth’s climate system. Current constructions of the global energy balance are based on 

38 the analysis of assimilated data, satellite estimates of global radiant energy and turbulent heat 

39 over oceans, and/or the hybrid approach of in-situ and satellite measurements [Da Silva et al., 

40 1994; Trenberth and Solomon, 1994; Rossow and Zhang, 1995; Yu et al., 1999; Trenberth and 

41 Stepaniak, 2004; Fasullo and Trenberth, 2007; Zhang et al., 2007; and references therein]. With 

42 these constructed atmospheric heat fluxes, atmospheric and oceanic poleward heat transports are 

43 estimated [e.g., Zhang and Rossow, 1997; Fasullo and Trenberth, 2007; Zhang et al., 2007]. 

44 Model assimilations can also provide global estimates of all atmospheric major energy and heat 

45 components. But significant errors associated with these estimates exist and can be as large as 

46 about 30 W/m" over large (1000 km) scales [ Trenberth and Solomon, 1994]. Some analysis 

47 techniques, especially the method of constraining the model analysis results with satellite top-of- 

48 atmosphere (TOA) radiation measurements and mass corrections within the assimilation models, 

49 are generally critical for reducing the uncertainties in global heat budgets [ Trenberth et al., 

50 2002], 

51 Satellite-estimated heat components of the global energy balance are mainly focused on 

52 the fluxes of TOA and surface radiative energy and air-sea turbulent heat [e.g., Wielicki et al, 

53 1996; Zhang and Rossow, 1997; Chou et al, 1997; Schulz et al., 1997]. Analysis of satellite data 

54 indicates that the mean differences among radiative flux data sets may be large enough that 

55 direct measurements of annual planetary energy imbalances are still unreliable. However, 

56 comparison of the interannual anomalies of the ocean heat content with satellite-derived 

57 planetary energy variations converted to accumulated ocean heat content (or equivalently 
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58 comparison of the anomalies of ocean heat storage converted from ocean heat content with the 

59 planetary energy imbalances) show excellent quantitative agreement [ Wong et al., 2006; Zhang 

60 et al., 2007]. Since both anomalies in and absolute values of the global energy budget are 

61 important for climate studies, quantitative knowledge about the global energy budget from more 

62 recent observationally-based data sets is needed. An earlier consistency study of blended 

63 satellite, in-situ and assimilation data for global annual mean atmospheric energy budget [Yu et 

64 al., 1999] found that the data sets available at that time resulted in an unbalanced atmospheric 

65 heat budget of 20 W/m , and the sign and magnitude of the systematic errors were consistent 

66 with the insufficient absorption of solar radiation within atmosphere debated at that time [e.g., 

67 Cess et al., 1995]. Although the systematic biases were generally much larger than TOA 

68 radiation uncertainties, these errors might be attributed to large spurious errors in the estimates of 

69 sea surface turbulent fluxes and to the combined effects of uncertainties in the radiation and 

70 turbulent flux calculations used in the study. 

71 Since there are significant improvements in both surface radiation and air- sea interaction 

72 flux estimates from satellite observations in last 5-10 years, this paper revisits the consistency 

73 issue of global annual atmospheric energy budget. The overarching goal is to evaluate the 

74 magnitude of the systematic biases within current satellite -based datasets and determine if the 

75 spurious errors are within the accuracies of current satellite retrievals of radiative and sea surface 

76 turbulent fluxes. The datasets are discussed in Section 2, and the results are shown in Section 3. 

77 Major conclusions are summarized in Section 4. 

78 

79 2. Data sets and analysis methodology 
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80 In this study, satellite observations are employed to estimate TOA radiative fluxes. For 

81 surface fluxes, satellite retrievals are used over oceans, and the combined results from satellite 

82 estimates of radiant energy and assimilation analyses of surface heat storage and the partition of 

83 latent and sensible heat (or the Bowen ratio) are used over land. Three global radiation datasets 

84 are used here: measurements from the Clouds and the Earth’s Radiant Energy System (CERES) 

85 mission [Wielicki et al. , 1996], the International Satellite Cloud Climatology Project Flux Data 

86 [ISCCP-FD, see Zhang et al., 2004], and the Global Energy and Water Cycle Experiment 

87 (GEWEX) Surface Radiation Budget (SRB) data [Stackhouse et al., 2001]. CERES directly 

88 measures TOA outgoing and incoming broadband longwave (LW) and shortwave (SW) radiation 

89 for the climate system. The other two radiation projects (ISCCP and SRB) calculate the TOA 

90 and surface radiation energy based on satellite observations of atmospheric temperature and 

91 humidity profiles, cloud optical properties and their spatial distributions, and the surface 

92 radiation properties such as skin temperature, emissivity and bidirectional reflection distribution 

93 functions. The random errors in the TOA monthly mean data at regional scales (-250 km) 

94 associated with these radiation data are reasonably small (-5 W/m 2 ; see the references listed 

95 above). The global monthly mean random errors are even smaller. The systematic errors in 

96 estimating the global annual mean energy budget can be as large as about 5 W/m" for the direct 

97 radiation measurements and within about 2 W/m" for ISCC-FD and SRB products. At the 

98 surface, the instantaneous errors in the radiative fluxes for the current ISCCP-FD and SRB 

99 products are as large as about 30 W/m". The regional monthly mean bias errors are significantly 

100 smaller, around 10 W/m" [ Zhang et al., 2004]. The system errors for global annual means could 

101 be even smaller due to potential cancellations of the bias errors for different climatological 

102 regimes. 
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103 The global turbulent heat fluxes from oceans to the atmosphere are based on the version 2 

104 and 3 products of the Goddard Satellite -based Surface Turbulent Fluxes (GSSTF) and Hamburg 

105 Ocean Atmosphere Parameters and fluxes from Satellite (HOAPS), respectively, and are 

106 estimated from satellite microwave sensors [Chou et al., 1997; Schulz et al., 1997]. The random 

107 error for instantaneous flux estimates is approximately 30 W/m , and that for monthly regional 

108 averages decreases to ~15 W/m . The systematic errors are much smaller and within about 7 

109 W/m'. Since there are no global land surface turbulent flux observations, the latent and sensible 

110 heat fluxes are calculated from a combination of the results from the Global Land Data 

111 Assimilation System (GLDAS) [ Rodlell et al., 2004] and the SRB radiation data. Because the 

112 temperature of regional land surfaces may vary from one month to another, there are small heat 

113 storage changes in the monthly time scale for a particular region. At the global annual mean 

1 14 scale, the land heat storage change [Huang, 2006] is much smaller than the systematic errors in 

115 the current datasets and the potential satellite-observed climate system energy imbalance. Our 

116 analysis confirms that the GLDAS yields negligible changes in the global annual mean heat 

117 storage. Also, the regional horizontal heat transport within land surfaces is much smaller than 

118 the storage change and can be ignored. Thus, this study uses surface SRB radiation and regional 

119 monthly heat storage from GLDAS as heat constraints for latent and sensible heat fluxes in each 

120 regional grid box (1.25°xl°). Furthermore, the monthly Bowen ratios in each grid box from 

121 GLDAS are used to partition the latent and sensible heat fluxes based on the heat constraints of 

122 SRB radiation and GLDAS storage fluxes. In this way, we have forced the land surface energy 

123 budget into balance at the global annual mean scales and essentially eliminate the spurious net 

124 flux errors over land. 
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125 Poleward of about 75°S, the surface is primarily covered by oceanic and continental ice 

126 sheets. There are few surface latent and sensible heat estimations from both satellites and 

127 GLDAS. Our satellite based estimates of global annual energy budget mainly cover the regions 

128 north of 75°S latitude. Because the turbulent fluxes are generally small south of 75°S, the 

129 sensible heat fluxes are assumed to be zero during cold seasons and the precipitation data from 

130 the Global Precipitation Climatology Project [GPCP; Adler et al., 2003] are used to fill the 

131 turbulent energy gap for these latitudes. Since the surfaces are very cold and there is only a 

132 small amount of moisture transported into the high latitudes, the latent heat estimated from 

133 precipitation and the assumed zero sensible heat fluxes from surface to atmosphere could 

134 overestimate the turbulent fluxes. On the other hand, due to GPCP underestimates of snowfall 

135 and drizzle, the overall errors in the estimates of the turbulent energy in the region may be 

136 reduced. Finally, all analyzed data are collected for the year 2000. There were no special 

137 climate events, such as significant El Nino, La Nina, or volcanic activities during this year. An 

138 analysis of that year’s satellite products represents the current status of satellite estimations of the 

139 global energy budget under normal climate conditions. Also, 2000 is the only year that satellite 

140 sea surface turbulent flux data from the GSSTF overlap with CERES radiation measurements. 

141 

142 3. Results 

143 Comparisons of the CERES, SRB and ISCCP TOA radiative fluxes reveal that the basic 

144 global patterns of annual mean TOA SW and LW fluxes, especially those for zonal averages, 

145 from all three data sets are very similar. The major differences are systematic biases among 

146 them, especially between CERES and the other two satellite calculations. As mentioned in the 

147 previous section, direct TOA radiation measurements yield a net radiation imbalnce of ~5.5 
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148 W/irf for the global annual mean, while SRB data result in a systematic imbalance of about 1.5 

149 W/m". Because this 5 W/m" imbalance has existed in the direct TOA radiation measurements for 

150 about 2 decades, it can be easily removed from interannual variation analysis, resulting in a 

151 much smaller (~0.5 W/m") residual systematic imbalance. In order to obtain a conservative 

152 annual energy budget and more realistic current satellite-based energy imbalance estimate, a 

153 somewhat larger bias in the SRB fluxes is considered here. Figure 1 shows zonal annual means 

154 of TOA (solid curve), surface (dotted curve), and atmosphere (dashed curve) net radiation 

155 estimates (note: hereafter all numbers in figures represent global mean values.) Integration of 

156 the TOA radiative fluxes from the poles to the equator represents the net meridional heat 

157 transports of the general circulation of the climate system. It can be seen from the TOA radiation 

158 plot that the climate system gains net energy only within - ±35° latitudes, and the middle 

159 latitudes have the maximum climate heat transports. The variation of zonal surface radiation 

160 basically follows the latitudinal pattern of TOA radiation except that the surface radiation is 

161 about 110 W/m" higher due to small differences in surface upwelling and downwelling LW 

162 radiation and to the dominant influence of solar radiation. The atmospheric net radiation, i.e., the 

163 difference between TOA and surface radiative fluxes is rather uniform, around -llOW/m" for 

164 most of latitudes. Within the atmosphere, SW absorption is minimal compared to LW emission 

165 and the LW radiation cooling into space dominates the atmospheric radiation budget. 

166 The annual zonal means of latent and sensible heat fluxes from the surface to the 

167 atmosphere estimated from GSSTF are shown in Figure 2. HOAPS produced results similar to 

168 those from GSSTF. Latent heat (solid curve) gradually decreases from more than 100 W/m" at 

169 low latitudes to nearly zero at poles. A clear relative minimum near the equator is caused by the 

170 weak winds of the intertropical convergence zone (ITCZ). Sensible heat fluxes (dashed curve) 
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171 are generally small compared to latent heat fluxes and range from about 0 to 25 W/m". The 

172 global annual averaged latent heat and sensible heat fluxes are 86 and 18 W/m , respectively. 

173 These latent heat fluxes are significantly greater (~ 1 1 W/m') than GPCP measured rainfall latent 

174 heat releases (dotted curve). Because there are basically no snowfall and drizzle estimates in the 

175 GPCP data set and significant uncertainties in both the rainfall and surface latent heat 

176 estimations, the two different estimates in the atmospheric latent heat are reasonably consistent. 

177 With full precipitation and surface latent flux retrievals, zonal moisture transports that currently 

178 have not been understood could be estimated. 

179 The annual zonal mean distribution of atmospheric total heat fluxes (Figure 3), the 

180 combined heating fluxes to the atmosphere from TO A and surface net radiation and surface 

181 latent and sensible heat, basically follows the latitudinal pattern of net radiation at TOA and 

182 surface except that a minimum exists at equator caused by the low surface turbulent heat fluxes 

183 at this region. Combining the strong atmospheric radiative cooling (112 W/m“) with the slightly 

184 weaker turbulent heat flux from surface to the atmosphere (104 W/m“), this analysis results in an 

185 estimated annual mean global atmospheric heat imbalance of about -8 W/m". Since the 

186 averaged atmospheric heat storage change in annual and global scales is negligible (considerably 

187 smaller than 1 W/m ), this global atmospheric heat imbalance is clearly a spurious error of the 

188 atmospheric heat budget. Similar to this atmospheric heat imbalance, the estimated global 

189 annual mean surface total heat imbalance is about 9.4 W/m". Although there has been some 

190 slight heating of the oceans and the earth’s climate system in recent years [ Wong et al., 2006], 

191 the relatively high value of 9.4 W/m in surface heating is largely the result of the various errors 

192 in the input data that caused a complementary bias in the atmospheric heat budget. When the 

193 systematic errors in turbulent (~7 W/m") and radiative (~10 W/m") heat fluxes are considered, 
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194 the systematic error (8~9 W/irT) in global total energy budget is not a surprise. Actually, this 

195 systematic error is less than half of what was estimated from the blended data of satellite, in-situ 

196 and assimilation in Yu el al. [1999]. Also, this spurious error is within the current understanding 

197 of the uncertainties in global radiation and turbulent flux estimates. Thus, there is no need to 

198 invoke the need for significantly more atmospheric absorption of solar radiation as mentioned by 

199 Yu et al. [1999] and as debated at about a decade ago. 

200 Global distributions of the oceanic annual mean surface heat budget are shown in Figure 

201 4. Positive values in the figure indicate that oceans gain heat from the atmosphere. Over land 

202 and at the annual time scale, there is almost no net heating due to the negligible heat storage and 

203 the forced balance among the radiative and latent and sensible heat fluxes, and the heat storage in 

204 this study, as mentioned before. Over oceans, regional net heating from the atmosphere is 

205 mostly used for horizontal heat transports with a relatively small part for vertical heat mixing. 

206 Since a portion of our estimates of the regional annual surface heat budgets, especially of those 

207 with small absolute numbers, is from bias errors in the regional estimations of radiative and 

2 

208 turbulent heat fluxes, the estimated annual budgets with an absolute value exceeding ~10 W/m 

209 could be significant for this analysis. For areas such as the ITCZ and those having strong ocean 

210 currents, heat horizontal transports dominate the estimated budgets. The equatorial area, 

211 particularly in the eastern parts of the ocean basins, is the major heat source of the oceans. It has 

212 a large net radiant energy gain, loses a comparatively small amount of turbulent heat, and has a 

213 surface heat budget as large as about 100 W/m 2 . The heat in the eastern ocean basins is 

214 generally moved to western basins by easterlies, then, transported to higher latitudes. Some of 

215 the surface heat to the ocean in these regions is also used for heating the upwelling cold water 


216 caused by the Ekman pumping. Both the Gulf Stream and Kuroshio Current play critical roles in 



217 latitudinal heat transports. They bring warm water from low latitudes to middle and high 

218 latitudes and release considerable latent heat into atmosphere. Combining turbulent cooling with 

219 radiative heating, we still find heat losses of more than 60 W/m" in these oceanic current regions. 

220 Large areas of the West Australia Current have cooling features similar to those of the Gulf 

221 Stream and Kuroshio Current except that the Australian current is much weaker. Oceans 

222 generally gain energy from the atmosphere over the annual time scale in tropical regions. 

223 Subtropical subsidence areas may have small annual heating budgets due to a combination of 

224 climate conditions of dry windy weather (i.e., large latent heat loss) and significant solar 

225 radiation. With rapidly decreasing in solar radiation with increasing latitude accompanied by 

226 smaller reductions in turbulent fluxes, the sea surface at higher latitudes releases heat into the 

227 atmosphere. It is because of the oceanic horizontal heat transport along with some vertical heat 

228 mixing, that the basic heat balance over sea surfaces is reached. The heat budget distribution in 

229 Figure 4 clearly shows major features of oceanic dynamics and the dominant mechanism of 

230 horizontal heat transports within oceans. 

231 

232 4. Summary 

233 This study uses the measurements taken in the year 2000 from multiple satellites to 

234 estimate global annual mean atmospheric heat budget. At the top-of-atmosphere, net radiative 

235 fluxes into the atmosphere obtained from both direct radiant energy measurements and radiation 

236 calculations using satellite-observed atmospheric profiles are close to zero. The global means of 

237 net radiative energy flux into the surface and surface latent heat flux into the atmosphere are 

238 about 113 and 86 W/m , respectively. The atmospheric and surface net heat budgets are about 

239 -8 ~ 9 W/m". These annual mean global heat imbalances in the atmosphere and at surface are in 
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240 the same order of magnitude as the uncertainties in the radiation and sea surface turbulent flux 

241 estimations and the likely systematic errors in the analyzed data. Although these spurious errors 

242 are significant for studies of annual mean global heat budget, they are clearly much smaller (less 

243 than half) than those estimated from blended data about decade ago [Yu el al., 1999]. 

244 Furthermore, the potentially strong additional absorption of solar radiation within the atmosphere 

245 as suggested by Yu el al. is not be required in the current analysis of the global energy budget 

246 due to much smaller spurious heat imbalances in the data compared to those used by Yu et al. . 

247 Progress in satellite surface radiation and oceanic turbulent heat flux estimations significantly 

248 reduces the bias errors in the observed global energy budgets of the climate system. 

249 Future work will be targeted on shrinking systematic errors in satellite estimates of 

250 surface radiative and turbulent heat fluxes. Removal of systematic heat budget errors would 

251 provide a great opportunity to use zonal annual means (such as those plotted in Figures 1 - 3) to 

252 estimate meridional heat transports of the earth’s climate system and separate the heat transports 

253 into atmospheric and oceanic components. Combining advanced precipitation measurements 

254 with surface latent heat estimations would also enable the estimation of atmospheric meridional 

255 moisture transports at an accuracy beyond that can be determined from the current, very limited 

256 measurements and observationally-based knowledge. 

257 

258 
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323 Figure captions 

324 Fig. 1. Annual zonal mean net radiation at TOA (solid), over surface (sfc; dotted) and within the 

325 atmosphere (dashed). Hereafter, the numbers for individual curves shown in the figure are 

326 their corresponding global annual means. 

327 Fig. 2. Annual zonal means of surface latent (solid) and sensible (dashed) heat fluxes. Also 

328 plotted is the latent heat (dotted) estimated from precipitation measurements. 

329 Fig. 3. Annual zonal means of atmospheric (solid) and surface (dashed) heat budgets. 

330 Fig. 4. Annual mean sea surface heat budget in W/m“. Positive values indicate that oceans gain 

331 heat from the atmosphere. 

332 

333 

334 
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335 Figures 


336 

337 Fig. 1 Annual zonal mean net radiation at TOA (solid), over surface (sfc; dotted) and within the 

338 atmosphere (dashed). Hereafter, the numbers for individual curves shown in the figure are 

339 their corresponding global annual means. 
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343 Fig. 2 Annual zonal means of surface latent (solid) and sensible (dashed) heat fluxes. Also 

344 plotted is the latent heat (dotted) estimated from precipitation measurements. 
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347 


Fig. 3 Annual zonal means of atmospheric (solid) and surface (dashed) heat budgets. 
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Fig. 4 Annual mean sea surface heat budget in W/m“. Positive values indicate that oceans gain 


heat from the atmosphere. 
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